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ABSTRACT 

We investigate the origins of galaxy morphology (defined by bulge-to-total K-band 
luminosity) in the ACDM cosmology using two galaxy formation models based on the 
Millennium simulation, one by Bower et al. (the Durham model) and the other by 
De Lucia & Blaizot (the MPA model). Both models have had considerable success in 
reproducing a number of observed properties of the local and high redshift universe, 
including star formation rates, the stellar mass function and the luminosity function 
out to z ~ 5. There are many similarities, but also fundamental disagreements in 
the predictions of the two models for galaxy morphology. For example, taking into 
account uncertainties in the available observational data, both produce a realistic 
morphological mix today, but its evolution is very different. A main cause of this and 
other differences is the treatment of disk instabilities which play a more prominent 
role in the Durham model. Our analysis confirms previous theoretical predictions that 
elliptical galaxies form most of their stars before the bulk of the galaxy is assembled. 
Spirals tend to have later 'assembly' times than ellipticals as a consequence of in-situ 
star formation. With the exception of the brightest ellipticals (stellar mass Af» ;> 
2.5 x 10 h Mq), we find that major mergers are not the primary mechanism by 
which most spheroids (ellipticals and spiral bulges) assemble their mass. In fact, the 
majority of ellipticals (and the overwhelming majority of spirals) never experience 
a major merger (above the resolution limit of our simulation.) Most ellipticals and 
spiral bulges acquire their stellar mass through minor mergers or disk instabilities. 
These conclusions are common to both the MPA and Durham models. The rotation 
properties of spheroids may help to constrain the importance of disk instabilities in 
these models. 

Key words: galaxies: formation - galaxies: evolution - galaxies: bulges. 



1 INTRODUCTION 

The morphological type of a galaxy describes more than 
merely its physical appearance; it reveals the key pro- 
cesses that have shaped it and continue to affect its evo- 
lution. H ubble's well e stabli shed classification scheme for 
galaxies (lHubblelll926l. Il93rj) h as survived with only mod- 
est modificatio n (|de Vaucouleursfl959l ; IVan den Berghfl960l ; 
ISandage] Il96ll ). precisely because the disks and spheroids 
on which it relies are fundamental products of the for- 
mation process. Understanding the origins of such struc- 
tures is therefore of central importance if we are to build 
up a coherent picture of galaxy formation and evolution. 
Modern observational facilities like the Hubble Space Tele- 
scope have made it possible to probe progressively higher 
redshifts and, together with large surveys like the Sloan 
Digital Sky Survey (SDSS), are beginning to reveal how 
galaxy morphology changes with time. A significant trend 
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in th e low redshift universe is the 'density-morphology' rela- 
tion i|Dresslerlll980l ; iPostman fc Gelled [l98l ; iDressler et "all 
1 1997T ) . which shows that denser regions contain proportion- 
ally more earl y-type galaxies. Thi s relation holds at least 
back to z ~ 1 I Pos tman et al . 2005). At higher redshifts, be- 
tween z ~ 0.5 and 2, there appears to be a rapid breakdown 
of the Hubble class ification system l|Van den Berghl 120021 : 
|Papovich~e t al. 2005) with the proportion of galaxies classed 
as irregulars or p eculiars rising steadily, par ticularly at faint 
magnitudes (e.g. jBrinchmann fc E llis 2000). Nonetheless, a 
substantial fractio n of giant spirals are obse rved at moder- 
ate redshifts (e.g.. iRavindranath et al1 l2004). Some authors 
have claimed that this population shows very little evolu- 
tion up to the present, implying that at least some Hubble 
sequence galaxies were in place at z ~ 1 (|Marleau fc Simardl 
1 19981 : IRavindranath eT~aTIl2004l ) . whilst others counter that 
many of the higher redshift sample would have become bulge 
dominated in the intervening time, so they are not di rectly 
comparable to their z = counterparts (Bell 2008). Be- 
fore 2 = 2, the majority of observed galaxies are either 
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highly irregular or centrally c ompact (|Giavalisco et al ] | 19961 ; 
iDickinsonS 2000; Daddi 2004). The difficulty in linking pop- 
ulations seen at different epochs has ensured that a clear 
evolutionary picture for morphology remains elusive. 

At present, direct simulation of the processes that de- 
termine morphology remains too computationally taxing 
to perform for large samples of galaxies (although signifi- 
ca nt progress is bein g made in this respect; see for exam- 
ple [Croft^etlaLlllQO^). Semi-analytic techniques constitute a 
convenient tool to capture the most important characteris- 
tics of these simulations, whilst also providing the statistics 
required to perform detailed analyses. The galaxy formation 
models considered in this paper are described in § [2] 

Our current understanding of galaxy formation tells us 
a little about how morphology develops. It is well known, 
for instance, that angular momentum acquired by a proto- 
galactic gas cloud through tidal torques operating during 
its expansion an d collapse will naturally c a use it to set - 
tle into a disk ijHovlel 1 19491 : IPeeblesl [i"9ol IWhitd 1 19841 ). 
Beyond this though, the combination of mechanisms re- 
quired tojHcplaffijnorphological change are poorly under- 
stood. iToomrd l|l977h was the first to put forward the idea 
that mergers could be responsible for the formation of el- 
liptical galaxies. Indeed, the importance of mergers in the 
hierarchical cold dark matter (C DM) model makes them a 
natural candidate for such a role ()Frenk et al ] 1 19851 ). Early 
CDM galaxy formation models stressed the possibility that, 
through a combination of mergers and gas accretion, a sin- 
gle galaxy could, at various stages in its lifetime, be iden- 
tified wi th several morpholog ical types on the Hubble se- 
quence l|White fc Frenkl 1991) , a view supported by gas- 



dynamical simulations Steinmetz fc Navarrd l|2002h . Numer- 
ous simulations have confirmed at least that major merg- 
ers, or very strong gravitational encounters, can lead to 
the complete destruction of a disk and the formation of a 
sphero i d with properties similar to e l liptic al galaxie s (e.g. , 
Whitd [l97Sl: iNegroponte fc White! 1 19831: Barnes! 1 19921; 
Barnes fc Hernauistlll996l:lBekki ^Tshioval[l997l: iNaab et all 
19991; ISpringel fc Herncmistl 120051 ; iNaab fc BurkertJ 120031 ; 



Cox et al.ll2006h . Even in virialised clusters where the rel- 



ative velocities of galaxies make strong gravitational inter- 
actions of any kind infrequent, mergers will still have af- 
fected groups collapsing to form the cluster, including those 
it accretes at later times. 

Other phenomena act to change morphology through 
their effects on gas content and star formation. For 
instance, strangulation (also called suffocation) occurs 
when a galaxy with a gaseous atmosphere falls into 
a larger stru cture, causing gas cooling to be heav- 
ily supressed (|Larson et all Il980l; IWhite fc Frenkl Il99ll ). 
In the formation mode ls of iBower et alj (120061 ) and 
|Pe Lucia fc BlaizotJ (|2007h considered here, strangulation 
is modelled by stripping all infalling (satellite) galaxies of 
their hot gas supplies, preventing further cooling. In com- 
bination with efficient feedback from supernovae, this can 
lead to the reddening of a satellite's stellar population on 
a relatively short timescale. Two further effects that are 
not i ncluded in these two models are ram p ressure strip- 
ping dGunn fc Gottj 1 19721 : lAbadi et all Il999h and harass- 
ment 1 Moore et al.l 19961 ). both of which are most impor- 
tant for cluster galaxies. The former involves the stripping 
of gas as galaxies move through the intracluster medium 



(ICM), while the latter results from a series of intense, high 
speed gravitational encounters that can cause bursts of star 
formation and lead to disk disruption. Recent additions to 
the Durham galaxy formation model analysed in this paper 
|Bower et al.1 120061) include a deta iled prescription for ram 
pressure stripping (|Font et al.| [2008). In this paper, we use an 
online database to examine two different galaxy formation 
models imp lemented on the ACDM M illennium Simulation, 
the first bv lDe Lucia fc Blaizot] il2007t) (here after the MPA 
model) and the second by Bower et al.l (|2006h (hereafter the 
Durham model). We look at the morphological content of 
each model as a function of redshift and track the evolution 
of the luminosity function of each type. We go on to deter- 
mine how a galaxy's formation epoch, and that of its host 
dark matter halo, affect its final morphology and investigate 
the role played by major mergers, minor mergers and disk 
instabilities. Our main goal is to investigate galaxies of all 
morphological types and to compare and contrast the pre- 
dictions of the two models, but we also carry out limited 
comparisons to observations. The rest of the paper is laid 
out as follows. In ijH we describe the N-body simulation and 
semi- analytic models that form the basis of this analysis. Sj3] 
discusses the methods we use to select and retrieve the data 
and compares some basic properties of morphology in the 
models with SDSS data. Finally, our results and conclusions 
are presented in Sj5]and i|6] 



2 THE N-BODY SIMULATION AND 
SEMI-ANALYTIC MODELS 

2.1 The Millennium simulation 

Both galaxy formation models investigated in this paper are 
based on the Millennium Simulation (MS) of the evolution 
of cold dark matter in a representative cosmological vol- 
ume. Completed by the Virgo Consortium in 2004, the MS 
remains the largest published N-body cosmological simula- 
tion, representing a significant step fo rward in terms of spa - 
tial, temporal and mass resolution (see lSpringel et all (|2005h 
for a detailed description). It uses 2160 3 (~ 10 10 ) particles 
to model a cubic region of space, 500/i -1 Mpc on a sideQ 
and tracks the evolution of structure between z — 127 and 
the present day. The initial configuration of the particles 
was constructed using a ACDM power spectr um consistent 
with cosmological parameter s from 2dFGR S ( Collcs s~et al.l 
1200 it) and first year WMAP (|Spergel et al.ll2003T l data. The 
parameters used were fl m = 0.25, fib = 0.045, h = 0.73, flA 
= 0.75, n — 1 and as = 0.93 at z = 0. Particle positions and 

velocities were recorded in 64 output snapshots. 

A friends-of-friends (FOF) algorithm (|Davis et al.l 

Il985h was used 'on the fly' to identify groups, linking ad- 
jacent particles separated by less than 0.2 times the mean 
interparticle separation. Groups with 20 or more particles 
(corresponding to a minimum halo mass of 1.7 x 10 10 /i _1 Mq) 



1 Here, h is defined by H = h X 100kms _1 Mpc _1 . 

2 Here, fi TO , fit an< i are the densities of all matter, baryons 
and dark energy respectively, in units of the critical density 
(pcrit = 3H 2 /8nG), n is the initial power spectrum slope and 
as is the rms overdensity predicted today by linear theory for a 
sphere of radius 8h~ ijvlpc. 
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were retained for further analysis. The FOF groups were 
further processed using a version of the SUBFIND algo- 
rithm (jSpringel et al.l I2001M that identifies gravitationally 
bound locally overdense regions, which to which we refer 
here as subhaloes. In the MPA model, a halo, by definition, 
includes all subhaloes within a FOF group. In the Durham 
model, on the other hand, a FOF group may be split into 
more than one halo if: i) a subhalo is outside twice the half- 
mass radius of the original halo, or ii) a subhalo was identi- 
fied in a previous snapshot as part of a different halo a nd has 
retained 75% or more of its mass |Harker et al.ll2006"h . Both 
of these conditions are designed to prevent two haloes that 
are temporarily joined by a tenuous particle bridge being 
linked together. The latter condition is applied on the basis 
that a subhalo loses significant mass from its outer layers 
when falling into a more massive counterpart, but substan- 
tially less if it is merely a close encounter. In both models, 
the descendant of a subhalo is identified by following its most 
tightly bound particles and the descendant of a halo is that 
which contains most of its most massive subhalo. With the 
merger trees of these bound structures thus defined, the two 
models employ different techniques to populate them with 
galaxies. For a detailed description of the two semi-analytic 
models analysed here, the reader is refe rred to the papers 
in which they were origin ally presented l|Bower et al.l 120061: 
|Pe Lucia fc Blaizoj 12007!) and references t herein , particu- 
larly ICroton et al.l i|2006l ) and ICole et~ai1 (|2000l ). The fol- 
lowing subsections summarise aspects that are particularly 
relevant to the discussion of morphology. 

2.2 Gas cooling, star formation and supernova 
feedback 

The treatment of radiative gas cooling in both models es- 
sentia lly follows the framework set out by IWhite fe Fre"nkl 
(1991). The amount of gas initially available to cool is a 
fixed fraction of the mass of the dark matter halc[f|. This gas 
shock-heats as it falls into the halo's potential well, attaining 
the virial temperature. As it cools and collapses, the angu- 
lar momentum that it acquired prior to turnaround causes 
it naturally to settle into a disk, which cools fastest in the 
centre where the density is highest. For low mass haloes, 
cooling is very rapid and the supply of cold gas is limited 
purely by the rate at which gas can free-fall onto the disk. 
For high mass haloes, a quasi-static hot atmosphere forms, 
from which gas can be accreted as it cools in the dense cen- 
tral regions. 

An important difference between the two models is in 
how they ca l culate the instantaneous rate of gas cooling. 
iBower et all {2006) allow gas to accrete onto the cold disk 
if it lies within a 'cooling radius', that at which the local 
cooling time is equal to the 'age' of the halo. They define 
the halo's age a s the time since it last do ubled in mass 
l|Cole et al.ll2000l ). |Pe Lucia fc Blaizod (|2007l ). on the other 
hand, calculate a rate based on the amount of gas that has 
a cooling time less than the halo dynamical time, which de- 
fines an alternative cooling radius. The Durham prescription 
typically results in higher cold gas masses than the MPA 

3 The reionisation of the gas at early times modifies this fraction 
somewhat for lower mass haloes in both models. 



prescription by a factor of ~ 2 at z < 5 in central galax- 
ies, rising to a factor of ~ 3.5 by z ~ 0. Stars a re assumed 
to fo rm from the col d disk gas, with Kennicutt ( Ke nnicutt! 
Il983h and Chabrier ([Chabrierl |20o3 ) initial mass functions 
in the Durham and MPA cases respectiv ely. In both models, 
the st ellar population synthesis model of lBruzual fc Charlotl 
|2003l ) is used to obtain stellar population properties. Star 
formation in the disk affects its development significantly. 
Stellar winds and supernovae increase the metallicity of the 
disk material from which further stars form and inject sub- 
stantial amounts of energy into their surroundings. This can 
drive gas and metals out of the disk, changing the compo- 
sition of the surrounding hot gas and altering its cooling 
time. 

2.3 Mergers and disk instabilities 

Two mechanisms incorporated into the Durham and MPA 
models are responsible for disrupting the stellar disks that 
would otherwise form as described above - galaxy merg- 
ers and disk instabilities. Galaxy mergers occur as a direct 
consequence of, but are distinct from, the mergers of the 
haloes that host them. As described in § 12.11 a halo that 
falls into a more massive system may survive for some time 
thereafter, as a gravitationally bound subhalo. Hence, the 
less massive (satellite) galaxy can be followed explicitly un- 
til tidal effects disrupt its subhalo sufficiently for it to drop 
below the 20 particle resolution limit. From this point on- 
wards, in the MPA model, the galaxy is associated with the 
most bound particle of the subhalo just before it became un- 
resolved. The satellite's orbit is assumed to decay through 
dynamical friction against the halo material until it merges 
with the more massive (central or primary) galaxy. Merger 
timescales are determined differently in the Durham model, 
with the satellite's orbit chosen at random fr om the cos- 
mological distribution derived by [Benson] l|2005l ). as soon as 
the halo merger has taken place and a dynamical friction 
timescale calculated accordingly. 

In the Durham model, the result of a merger is dictated 
by the relative masses of the merging galaxies and the gas 
content of the primary. 'Major' mergers (M aa t/M pr i ^ 0.3), 
completely disrupt any existing disk, producing a spheroidal 
remnant that contains all the stars from its progenitor galax- 
ies. Any gas present forms stars in a burst, which are also 
added to the new spheroid, with some fraction of gas be- 
ing returned to the hot halo through supernovae feedback. 
Minor mergers (M 3a t/M pr i <0.3) leave the structure of the 
primary intact, simply adding the stars from the satellite to 
the bulge of the primary and its gas to the disk, but may 
also initiate a burst if M sa t/M pr i ^0.1 and the primary has 
sufficient gas in its disk (M g as/Mdisk >0.1). The fraction of 
the available gas that is turned into stars is dictated by star 
formation and feedback rules, applied over the timescale of 
the burst. 

The MPA model assumes starbursts to occur in all 
mergers. The mass fraction of gas available to take part 
is a power-law function of the satellite-central mass ratio 
(m s at I'm central)- Other than this, the mechanics of mergers 
are largely similar to the Durham model. Major mergers are 
defined in the same fashion and also result in the forma- 
tion of a spheroid; minor mergers add stars to the primary's 
bulge and gas to its disk. 
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In galaxies with strongly self-gravitating disks, that is, 
where the mass of the disk itself dominates the gravitational 
potential, the system is unstable to pe rturbations from 
nearby satellites or d ark matter clumps (|Efstathiou et all 
Il982l ; IMo et all 119981 ). Both models regard a disk to have 
become unstable when 

V max /(GM Msk /r dlsk ) 1/2 < 1, (1) 

where Mdisk and r^isk are th e disk's total mass and radius. 
In the original formulation bv lEfstathiou et all (|l982h . V max 
was the maximum of the rotation curve, but this is approxi- 
mated by the halo virial velocity in the MPA model and the 
disk's velocity at its half mass radius in the Durham model. 
Although the prescriptions will typically agree on whether 
a particular disk is unstable or not, they deal with the sit- 
uation somewhat differently. In the MPA model, there is 
effectively a partial collapse, intended to model the forma- 
tion of a bar. Mass is moved from the disk into the spheroid 
until the stability of the system, defined by eqn. [1] is re- 
established. In the Durham model, instabilities result in the 
complete collapse of the disk into a spheroid, the size of 
which is dictated by the rotational energy of the disk just 
prior to collapse. As with a major merger, a starburst is 
induced, such that the resulting spheroid contains the orig- 
inal disk's stars, plus those formed in the burst. This more 
catastrophic outcome is assumed to result from orbital res- 
onances and stellar scattering in the barred system causing 
it to collapse entirely. 



where a quasi-static hot halo has formed, the energy gener- 
ated through accretion is assumed to couple efficiently with 
the hot halo gas. If the energy output of the BH (some con- 
stant fraction of its Eddington luminosity) exceeds the rate 
at which the gas can radiate away energy, then no further 
gas is allowed to cool. Hence, as soon as a large cooling flow 
builds up, the feedback becomes sufficient to cut it off. This 
behaviour plays a key role in reproducing the observed K 
and B-band luminosity functions out to z ~ 5. 

The MPA implementation also assumes that gas in- 
flows associated with mergers lead to increased BH accretion 
rates. Consistent with their starburst model, such inflows are 
presumed to occur in all mergers, with the accreted fraction 
of the total gas dependent on the satellite-central mass ra- 
tio. A simple phenomenological model is used to describe 
accretion from quasi-hydrostatic cooling flows; the mass ac- 
cretion rate (m) depends on the mass of the black hole, the 
hot gas fraction and the virial velocity of the halo. The lu- 
minosity of the BH at any given time is simply r/rhc 2 , where 
rj ~ 0.1 is the efficiency of mass-energy conversion typically 
expected close to the event-horizon. Based on this energy 
output, an adjusted cooling rate is calculated, such that it 
declines smoothly as the amount of mass accreted rises. As 
in the Durham model, BH-BH mergers occur and are mod- 
elled in a similar fashion. However, no additional accretion 
is considered during a disk collapse. 



2.4 Black holes and AGN feedback 

Prescriptions for supermassive black hole (BH) growth and 
feedback from active galactic nuclei (AGN) are relatively 
recent additions to both models, motivated by the grow- 
ing body of evidence linking propertie s of a galaxy's bulge 
to the mass of its central BH (e.g ., ?Fcrrarcsc fc Merritj 
20001; iKormendv fc Gebhardtll200ll ; iMarconi fc Huntll2003l ; 
Haring-Neumaver et alj|2006l ). Energy produced through ac- 
cretion of material onto the BH is assumed to restrict gas 
cooling, particularly in higher mass haloes, impacting on 
galaxy luminosity, colour and morphology. Different imple- 
mentations of AGN feedback could, for instance, dictate 
whether or not an elliptical galaxy is able to regrow a stellar 
disk after a major merger. 

The prescription used in the Durham m odel for AGN 
feedba ck and BH growth is described in full in lMalbon et al.l 
(|2007h . In the Durham model, BH growth is assumed to 
take place through four main channels, i) Major merger- 
driven accretion - Tidal forces act to drive disk gas into 
central regions, fuelling the BH. A constant fraction of the 
accreted gas, Jbh, (tuned to match the amplitude of the 
Mbh — Msuige relation) is added to the BH. ii) Instability- 
driven accretion - A fraction, Jbh, of the gas in the collaps- 
ing disk is added to the BH. iii) BH-BH mergers - When two 
galaxies merge, their BHs are assumed to merge. Mass loss 
due to the radiation of gravitational waves is neglected and 
the mass of the new BH is the sum of the progenitor masses 
plus any gas accreted, iv) Accretion from cooling flows - In 
sufficiently massive haloes, where quasi-hydrostatic cooling 
is taking place, the BH can accrete mass from the cool- 
ing flow. This last mechanism is closely associated with the 
implementation of AGN feedback in the model. In haloes 



3 THE DATA 

3.1 The Millennium Run database 

All of the data analysed here comes from the online Mil- 
lennium Run databasfl developed by the German As- 
trophysical Virtual Observatory (GAVO), which can be 
accesse d through the use of structured query language 
(SQL) (|L cmson fc The Virgo Consortium! |2006h and con- 
tains galaxy and halo data from the Durham and MPA 
models. Various properties are associated with each halo, for 
instance, mass, redshift, position, velocity, number of sub- 
haloes, etc. Galaxies additionally carry many directly and 
indirectly observable properties, such as observer and rest 
frame magnitudes in several passbands, stellar mass and ra- 
dial extent of morphological components. The database has 
been constructed to make the ret rieval of an object's merger 
tree as efficient as possible (see iLemson fc Springell (2006) 
for a detailed description). This has been achieved by as- 
signing unique IDs in a 'depth-first' manner, such that the 
progenitors of a given object have IDs lying between that of 
the object itself and an index it carries called the lastprogen- 
itorid. Additional indices identify each object's immediate 
descendant and largest progenitor. An additional index in 
the Durham tables identifies every galaxy's main branch. 



4 Hosted at http://www.g-vo.org/MyMillennium with a mirror 
site at http://galaxy-catalogue.dur. ac.uk:8080/MyMillennium 
(Both require registration for full access). 
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3.2 Definition of morphology and construction of 
galaxy samples 

To describe a galaxy's morphology we use the bulge-total 
(B-T) ratio of absolute, K-band, rest frame luminosity 
(R — L k, bulge / Lk, total)- Although most photometric mea- 
surements are made in the B-band, we chose the K-band be- 
cause it reflects th e stellar mass quite closely, even at moder- 
ate r edshifts (e.g. iKauffmann fc Chariot! 1 19981 : lLacev et al.l 
2008). This choice has the further advantage that the pre- 
dictions are then relatively insensitive to uncertain details of 
the current star formation and reddening corrections in the 
model. Galaxies are divided into three broad morphological 
types: "spirals" (R < 0.4), "SOs" (0.4 < R s: 0.6) and "ellip- 
ticals" (R > 0.6). This classificati on is, to some exte nt, ar- 
bitrary, but at least in the B-band [Tran et ail (|200lT ) found 
that galaxies classified as late Hubble types in the NASA 
Extragalactic Database are well by a rest-frame B-band B- 
T ratio of less than 0.4, consistent with our definition of 
spirals. 

In addition, we split our two z — galaxy populations 
by luminosity, with the division at Mk — 51ogh = —22.17, 
i.e. one magnitude fainter than the characteristic lum inos- 
ity in the K-band l|Cole et alj|200ll : ISmith et al.ll2008l . e.g.). 
Throughout this work, we refer to these as the bright and 
faint populations. Where the r-band is also shown, we di- 
vide the populations at M r — 51ogh = —19.83 jBlanton et al.l 
1200 ll ) and define morphological c lasses according to the same 
B-T ratios as in theK-band. iDe Lucia fc Blaizod (|2007l ) 
note that, in order to assign a morphology to a galaxy 
with confidence, its merger history must be well resolved. 
They determine that this condition imposes a lower limit of 
4 x 1O 9 /i -1 M0 in stellar mass, which we apply consistently 
to the bright and faint populations in both models. This cut 
has virtually no impact on the bright population but it re- 
duces the numbers in the faint population by ~ 92 — 94%, 
leaving us with 1,298,118 bright and 1,889,131 faint galax- 
ies in the MPA model and 1,280,154 bright, 1,598,908 faint 
galaxies in the Durham model. 



4 OVERVIEW OF MORPHOLOGICAL 
CHARACTERISTICS 

4.1 Evolution of the morphological mix 

With our galaxy samples thus defined, we now consider some 
basic properties of the two models. Firstly, we examine their 
morphological content, that is, the relative fractions of each 
morphology at a given redshift. Fig. [T] tracks the fraction 
of the total bright galaxy populations in the K-band (left 
panel) and r-band (right panel) that are spirals, SOs and 
ellipticals in each model. The r-band is includ ed here in or- 
der to compare with the observational data of iBenson et al.l 
l|2007l ). 

Clearly, in both bands, the two models differ consider- 
ably, particularly at high redshifts. The Durham model pre- 
dicts a substantial population of ellipticals at early times, 
with steady evolution to a disk-dominated phase after z ~ 
3 — 4. In contrast, the MPA model shows more modest evo- 
lution, with spirals remaining the most prevalent structures 
throughout. Even so, the two prescriptions result in fairly 
similar present day morphological compositions (see table [1] 



Table 1. Percentage contribution of each morphology to the 
present day populations of the MPA and Durham models. 
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for K-band fractions). Fractions in the Durham model are 
consistent with the ranges defined by the observational data, 
but the MPA model seems to produce too many spirals and 
not enough SOs. We will explain the origin of these ranges 
in our discussion of Fig. [2] The difference between bands for 
both models is a small offset toward higher B/T values in 
the r-band, which is a consequence of an increased sensitiv- 
ity to dust obscuration compared to the K-band, which acts 
to dim disk light. 

Considering instead the faint galaxy population (Mk — 
51ogh > —22.17, not shown here), we find a very similar 
result. The plot looks almost identical for MPA galaxies and 
similar for Durham galaxies, though there is an offset of 
about 10 — 15% in favour of more spirals at the expense of 
SOs and ellipticals. 

The disagreement at high redshift evident in Fig.[T]may 
be explained by the contrasting treatment of disk instabil- 
ities in the two models. As outlined in § 12.31 instabilities 
trigger a total collapse of the galactic disk in the Durham 
model, but only a partial "buckling" in the MPA model. As 
some of the results in § [5] demonstrate, this distinction ap- 
pears to be the root cause of several differences in galaxy 
morphology in the two models. 

Several attempts have been made to quantify morphol- 
ogy through bulge-total measurements at low redshift (e.g. , 
i Tasca fc Whitell2005l ; IBenson et al.ll2007l; iDriver et al.ll2007l : 
iGado tti 2008). In partic ular, Benson et al.l l|2007r ) showed 
that the Durham model (|Bower et al.ll200q) reproduces the 
luminosity function of disks and spheroidqj remarkably ac- 
curately. Fig.[2]uses their SDSS data to compare the fraction 
of spirals, SOs and ellipticals in the local universe (z < 0.12) 
as a function of (AB, rest-frame) r-band magnitude, with 
fractions from the two models at z = 0. For this plot, mor- 
phologies in the models are defined using rest-frame r-band 
light for consistency with the data. 

Benson et al. determined B/T ratios for their sample of 
~ 9000 SDSS galaxies by fitting two component light profiles 
(bulge + disk) to each image. Analysis of their initial dataset 
showed that the distribution of disk inclination angles ap- 
peared to be biased, leading them to conclude that their code 
was fitting face-on disks where the bulge varied significantly 
from the anticipated De Vaucouleurs profile. They applied a 
statistical correction which attempted to make the distribu- 
tion of inclination angles uniform, setting B/T to 1 for those 
galaxies left over after the correction. As Benson et al. point 
out, this correction will likely overestimate the prevalence 
of ellipticals and underestimate the number of spirals, but 
the magnitude of the uncertainty is very difficult to quan- 
tify. We have therefore chosen to show the data before and 
after the applied correction to illustrate the range of un- 



5 Spheroids were taken to include both the bulges of spirals and 
elliptical galaxies. 
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Figure 1. The fraction of bright galaxies classified as spirals (blue), SOs (green) and ellipticals (red) in the Durham (solid lines) and 
MPA (dashed lines) models as a function of redshift. In the left panel, morphology is defined using rest-frame K-band light and in the 
right panel using rest-frame r-ba nd light. The bars o n the right of the r-band panel represent the possible range in each fraction at 
2 = 0, calculated using the data of lBenson et alj l|2007f ). We explain the origin of these observed ranges in our discussion of Fig. [2] Disks 
are the dominant structures in the MPA model at all times, whereas the Durham model produces a transition from a bulge-dominated 
phase before z ~ 3 — 4 to a disk-dominated phase therefafter. This strikingly different behaviour appears to arise from the contrasting 
treatment of disk instabilities in the two models. Despite the differences, the two prescriptions arrive at fairly similar fractions at z = 0. 




M r -f>logh M r -5logh M r -5logh 



(a) (b) (c) 

Figure 2. The fraction of galaxies at z = that are spirals (a), SOs (b) and ellipticals (c), a s a function of rest- frame r-band magnitude 
in the MPA model (thick dashed line), Durham model (thick solid line) and in the SDSS data ([Benson et al,ll2007ri (connected points). The 
two sets of points represent the data before (triangles) and after (squares) a correction for biases introduced in the fitting process. This 
correction is explained in the text. The upper and lower magnitude limits of the plot are dictated by the data. Both models reproduce 
the general trend and typically lie between the uncorrected and corrected data, although the MPA model appears to have too many 
spirals and not enough SOs over a range of magnitudes. 



certainty for each fraction. The ranges at z — in Fig. \T\ 
are obtained by integrating the corrected and uncorrected 
fractions f or each type, multip lied by the r-band luminosity 
function (|Blanton et alj|200ll ). 

With this in mind, we see in Fig. [5] that both mod- 
els reproduce the general trend in the data - ellipticals are 
more prevalent at brighter magnitudes, and typically fall be- 
tween the corrected and uncorrected datasets. The fraction 
of spirals in the MPA model exceeds the uncorrected data at 
intermediate magnitudes and the fraction of SOs falls below 



the corrected data over the same range, but by an amount 
less than size of the Poisson errors on each datasets (not 
shown) . Consequently, neither model is obviously discrepant 
with the data, although clearly smaller uncertainties would 
be desirable in order to place tighter constraints. 
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4.2 The morphology-dependent luminosity 
function 

The luminosity function provides an important census of 
any galaxy population. Fig. [3] shows the rest-frame B and 
K-band luminosity functions for each morphological type 
at z = 0, the solid lines representing the Durham model 
and the dashed lines the MPA model. A galaxy's K-band 
luminosit y can be used as a reasonable proxy for its stel- 
lar mass dKauffmann fc Charlot|[l99Sl ; iBrinchmann fc Ellis! 
l200d ; iBell fc de Jonell200lf ). whilst the B-band is biased to- 
ward younger stellar populations and hence is indicative 
of recent star formation. On both plots, the black vertical 
dashed line marks the magnitude faintward of which the im- 
posed mass cut at 4 x 10 9 h -1 Mq depletes galaxy densities 
by > 5%. 

Both models produce luminosity functions which, qual- 
itatively at least, share some features with comparable ob- 
servational data. For instance, the most massive (brightest 
K-band) galaxies are ellipticals a nd spirals have a some what 
fainter characteristic mag nitude l|Nakamura et al.ll2003l '). We 
also note that spiral galaxies in the MPA model are typically 
brighter in the B-band than those in the Durham model, 
suggesting that they have significantly more ongoing star 
formation or less dust obscuration. In both models, ellipti- 
cals are the most common morphology at high (M < M*) 
K-band luminosities, but at the bright end of the B-band 
LF, the MPA model has at least an order of magnitude more 
more spirals than the Durham model. 

Fig. |3] shows the redshift evolution of the K and B- 
band luminosity functions for each morphology as ratios to 
the z = functions in Fig. [3] The reader should note that 
the bright-end extent of these curves is limited by that of 
the z = data and the faint-end is cut off where the im- 
posed mass cut reduces number densities by > 5%. In the 
K-band, there is very little evolution in the Durham model 
between z ~ 2 and z ~ 0.5 for all spirals and SOs and for 
all but the brightest ellipticals, suggesting that a substantial 
proportion of their mass was already in place at that time. 
At the bright end, there is already a sizeable population of 
ellipticals at z — 5 in both models. The early dominance of 
ellipticals in the Durham model seen in Fig.[T]comes mostly 
from galaxies excluded from this plot by the faint-end cut 
off. The number density of bright ellipticals in the Durham 
model declines slightly from z ~ 2 to z ~ 0.5, a result both 
of mergers between ellipticals and their transformation into 
SOs or spirals for those able to grow a disk. The trend is 
in the opposite sense for the MPA model, indicating that 
the extra spirals already in place at z — 5 compared to the 
Durham model make merging a process that typically in- 
creases elliptical number density. Evolution of the faint-end 
K-band is more significant in the MPA model, particularly 
for ellipticals. This may reflect weaker supernovae feedback 
allowing more star formation in satellite galaxies. 

The Durham B-band functions are consistent with a 
peak in star formation somewhere around z = 2, for spirals 
and SOs, and somewhat earlier for ellipticals. In the MPA 
model, there is a steady decline in B-band light for all types 
since z ~ 5, suggestive of consistently earlier star formation 
than in the Durham model. The history of star formation 
in the different types is examined more closely in the next 
section. 



5 RESULTS AND DISCUSSION 

In the previous section, we saw that there is some disagree- 
ment between the Durham and MPA galaxy formation mod- 
els in their predictions for the evolution of the relative abun- 
dances of galaxies of different morphological type. In this 
section, we explore the reasons for these differences by inves- 
tigating the formation paths of galaxies in the two models. 
Firstly, we track the build up of a galaxy's dark halo and its 
stars with the goal of identifying how its final morphology is 
determined. We then examine in more detail how spheroids 
form in each model and the roles played by mergers and disk 
instabilities. 



5.1 Halo formation 

As described in § 12.31 galaxy mergers follow the mergers of 
their host haloes after a time delay due to the gradual decay 
of the satellite's orbit through dynamical friction. One might 
expect, then, that the manner of a halo's formation might 
play a significant role in determing the merger history of its 
galaxy and hence its final morphology. 

We define the formation redshift of a halo, Zhform, as 
that at which half of its final mass was contained in a sin- 
gle progenitor. Haloes in the database are linked, through 
their unique IDs, to the galaxies that they host, making it 
straightforward to relate the time of formation of a halo 
to the morphology of its galaxies. We consider only 'FOF 
haloes' since they can be associated with galaxies in both 
models, as opposed to the modified Durham haloes which 
cannot always be linked to MPA galaxies. This allows us 
to compare the morphologies of galaxies in the two mod- 
els for a common halo population. We consider only central 
galaxies and limit our sample to haloes in the mass range 
0.5 — 2 x 10 12 h~ 1 Mo, a factor of two either side of the mass 
of a 'Milky Way type' halo. 

Fig. [5] shows the distribution of zh form for bright spi- 
rals, SOs and ellipticals in the Durham (solid line) and MPA 
(dashed line) models. Relatively few galaxies in the faint 
sample exist in haloes within this mass range, making the 
equivalent plot extremely noisy. 

Given that major galaxy mergers disrupt disk struc- 
tures, one might naively have expected later forming haloes, 
which merged more recently, preferentially to host ellipti- 
cals. In the Durham model, this effect is marginal, but can 
be seen in the median formation redshifts for haloes hosting 
each type. Halos with elliptical central galaxies have a me- 
dian formation redshift of 1.17, compared to 1.28 for those 
with spiral or sO centrals. In the MPA model, haloes host- 
ing ellipticals have a median formation redshift of 1.28, in 
comparison to 1.38 for those hosting SOs and 1.17 for those 
hosting spirals. 

The spike in haloes forming near z — 0, seen in both 
models, appears to arise from the difficulty in tracking ob- 
jects between snapshots. When a halo passes through a more 
massive counterpart, it may appear to have merged with it, 
but, if it emerges on the other side, it is labelled as a newly 
formed object. In this way, some merger tree branches are 
broken. If the break occurs after the real formation time 
and on the main branch, the final time halo will appear to 
have formed at that epoch. In most cases, the less massive 
halo falls back into its counterpart and merges before the 



8 0. H. Parry et al. 




last snapshot; thus it is no longer considered a main branch 
progenitor. Hence, the spike at z = is a boundary effect 
- it identifies those cases where insufficient time has passed 
for a subsequent merger to erase evidence of a fly-through. 

The above explanation is supported by the fact that the 
late-forming haloes are invariably found close to the virial 
radii of larger haloes. Furthermore, we find that the low 
redshift spike disappears if we consider instead the formation 
of the Durham haloes, for which considerable effort has been 
made to r epair broken merg er trees and eliminate such 'false 
mergers' l|Hellv et alj|2003f ). 



5.2 Galaxy assembly and formation 

Having established the extent to which galaxy morphology 
depends on the epoch of halo formation, we now consider 
the formation of galaxies directly. We define two different 
'formation epochs': i) the stellar assembly redshift, z a> when 
half of the stars in a present-day galaxy were first assembled 
in a single progenitor, and ii) the stellar formation redshift, 
Zf, when half of the stars in a present-day galaxy had formed 
in any of its progenitors. The three panels in Fig. [5] compare 
the distributions of these two redshifts for spirals, SOs and 
ellipticals. 

The two models show the same general trend for both 
formation and assembly to occur earliest in elliptical galaxies 
and latest in spirals, with SOs intermediate between the two. 
Given that early-type galaxies are typically observed to have 
older stellar populations, it is reassuring that both models 
predict that they form most of their stars first. Less intuitive, 
perhaps, is the fact that ellipticals also assemble their stellar 
mass before either SOs or spirals. Late assembly, however, 
does not imply late merging. Steady, in-place star formation 
in a disk naturally results in a relatively narrow distribution 
of assembly times, which are typically later than for elliptical 
and SO galaxies where mergers play a more important role. 
The fact that the formation and assembly epochs of spirals 
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Figure 5. The distribution of formation redshifts for haloes host- 
ing galaxies of each morphological type in the Durham (solid) and 
MPA (dashed) models. Only haloes in a 'Milky Way type' mass 
range (0.5 — 2 X 10 12 h~ 1 Mq) are considered here for the reasons 
outlined in the text. The median formation redshifts for spiral, 
SO and elliptical hosts are 1.17, 1.38 and 1.28 respectively in the 
MPA model, and 1.28, 1.28 and 1.17 in the Durham Model. All 
curves are normalised such that the area under them is unity. 



(and also SOs in the Durham model) are generally coincident 
suggests that this is indeed the case: when the formation 
threshold is reached, the stars are already in one object. 
In the Durham model, disks typically form at z ~ 0.55, 
although a significant fraction form at z > 1. In the MPA 
model the typical formation redshift of disks is z ~ 0.9. 

In both models, we effectively see two characteristic 
formation times for ellipticals and SOs, a primary peak at 
z ~ 1.5 — 2 and a secondary peak at z ~ 0.55. The latter co- 
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Figure 4. The evolution of the K-band (left column) and B-band (right column) rest-frame luminosity functions for spirals (top row), 
SOs (middle row) and ellipticals (bottom row), plotted as ratios to the z = luminosity functions in Fig. [3] The curves are cut off at a 
magnitude faintward of which the mass cut at 4 X 10 9 fr — 1 Mq reduces number densities by > 5%. At the bright end, they are cut off at 
the limit of the equivalent z = function. 
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Figure 6. The distributions of the stellar assembly and stellar formation redshifts for bright (Mr — 51ogh < —22.17) spirals (a), SOs 
(b) and ellipticals (c) in the MPA (dashed line) and Durham (solid line) models. All curves are normalised such that the area under 
them is unity. 



incides with the typical epoch of spiral galaxy formation and 
is more pronouced in the Durham model. This secondary 
peak is the result of disk instability: these ellipticals and 
SOs are, in effect, failed spirals, having had very similar star 
formation histories up to z ~ 0.55, before becoming unstable 
and collapsing into a spheroid with the attendant conversion 
of residual gas into stars. Some of these spheroids may par- 
tially regrow a stellar disk and end up as SOs or perhaps even 
as spirals, but others remain as ellipticals up to the present 
day, either because the instability occurred very late on, or 
because their gas cooling rate was restricted by feedback or 
lack of gas. Since disk instabilities have a much gentler ef- 
fect in the MPA model, their influence on the distribution 
of stellar formation times of ellipticals and SOs is minor. 

Minor and major mergers play an important role in 
building up most of the ellipticals in the MPA model and 
those whose stars form in the higher redshift peak in the 
Durham model, as evidenced by the fact the stellar forma- 
tion times are earlier than the stellar assembly times. This 
asynchronism between stellar formation and assembly in 
C DM-based models o f elliptical galaxies was first highlighted 
by iKauffmannl (1 199(f) and has recen t ly bee n investigated in 
the MPA model by |Pe Lucia et~aT1 (|2006l ). Our results for 
MPA ellipticals are in agreement with theirs. The distribu- 
tions of z a and 2/ for faint galaxies (Mk — 51ogh < —22.17) 
in both models are similar to those for bright galaxies but 
shifted slightly to higher redshift and with a more signifi- 
cant high redshift tail. Also, the lower redshift peak in the 
Durham model, corresponding to galaxies forming through 
disk instabilities, is much less pronounced for faint objects. 

5.3 Major mergers, minor mergers and disk 
instabilities 

The differences between the Durham and MPA models un- 
covered so far appear to stem from the roles played by the 
main spheroid-forming mechanisms. In this section, we com- 
pare the relative contributions of three processes that con- 
tribute to the build-up of spheroids: major mergers, minor 
mergers and disk instabilities. We remind the reader that in 



both models the distinction between major and minor merg- 
ers is made on the basis of the relative baryonic masses of 
the satellite and primary, with the boundary corresponding 
to M S at/Mpri = 0.3. In both models, major mergers are as- 
sumed to disrupt completely any pre-existing galactic disks. 

The incidence of major mergers and the redshift when 
the main branch progenitor of a present-day galaxy was 
last involved in a major merger, z mm , can be readily ob- 
tained by tracing back the galaxy's merger tree. The anal- 
ysis below includes galaxies above our resolution limits of 
1.7 x 10 10 h^ 1 Mq for the dark matter mass of the halo and 
4 x 10 9 7i _1 Mq for the stellar mass of the galaxy. The statis- 
tics of major mergers are surprising and revealing. Of the 
total bright M K - 51ogh < -22.17) population, only 49% of 
ellipticals, 2% of SOs and 2% of spirals in the MPA model, 
or 41%, 2.5% and 1.5% respectively in the Durham model, 
undergo a main branch major merger in their entire forma- 
tion history. In both models, this fraction is largely inde- 
pendent of total stellar mass for spirals and SOs, but not so 
for ellipticals. Combining the bright and faint populations 
together, in the Durham model, 100% of the highest mass 
(1 - 4 x 10 12 /i _1 M Q ) ellipticals have major mergers, but 
this fraction falls to < 3% around the lower mass limit of 
4 x 10 9 h _:L M©. In the MPA model, only ~ 15% of ellipticals 
with mass 1.6 x 10 10 /i _1 Mq have a major merger, but the 
proportion rises sharply either side of this mass, to around 
~ 80% at the low mass end, and close to 100% for masses 
> 2.5 x 10 n /i _1 M©. In both models, more than half of the 
total elliptical population do not undergo a main branch ma- 
jor merger at any point. These results are consi stent with 
the conclusions reached bv lDe Lucia et al.l (|2006h for MPA 
ellipticals. 

The distribution of z mm for each morphological type in 
the bright population is displayed in Fig. [7] There is rea- 
sonable agreement between the two models for spirals and 
ellipticals, but not for the SOs, which appear to trace the 
elliptical distribution in the MPA case, but the spiral dis- 
tribution in the Durham case. As one might expect, ellipti- 
cals typically have had the most recent major mergers, but 
nonetheless, of the small fraction of spirals and SOs that 
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Figure 7. The distribution of the last (main branch) major 
merger experienced by spirals (blue), SOs (green) and ellipticals 
(red) in the Durham (solid line) and MPA (dashed line) bright 
{M-K — 51ogh < —22.17) galaxy populations. Each curve is nor- 
malised such that the area under it is unity. The median values 
of Zmm are 2.62, 0.83, and 0.41 for spirals, SOs and ellipticals 
respectively in the MPA model, and 2.07, 1.63 and 0.51 in the 
Durham model. 



have undergone major mergers at all, some experience them 
at very late times (z < 0.1) and yet still manage to recover a 
stellar disk. The faint population shows similar behaviour to 
that illustrated in Fig. [7] Given that major mergers are this 
infrequent, the alternative mechanisms of minor mergers and 
disk instabilities must be important in spheroid formation. 
Applying a similar analysis to instabilities reveals that vir- 
tually all ellipticals have experienced such an event at some 
point in their formation history. However, the fraction of the 
total galaxy population that is unstable at a particular red- 
shift differs substantially between the two models. In any 
one snapshot, approximately 1% of MPA galaxies are un- 
stable but, in the Durham model, the figure is much higher: 
~ 16% at 2 = 6, falling to ~ 2% at z = 2. This difference 
accounts for the higher proportion of ellipticals seen at early 
times in the Durham model, as illustrated in Fig. [T] Early 
forming disks in the Durham model are inherently less sta- 
ble than those in the MPA model, largely as a consequence 
of their higher cold gas content. 

To quantify the relative impact of instabilities, major 
mergers and minor mergers on galaxy morphology, we de- 
termine the fraction of stellar mass in present-day spheroids 
that was incorporated into the spheroid by each process. 
The complexity of the SQL query needed for this calcula- 
tion made it necessary to use the 'milli-millennium' datasets, 
which were generated using the same two semi-analytic mod- 
els on halo merger trees constructed from a smaller version 
of the Millennium simulation (1/512 of the volume). Even 
though the milli-millennium simulation has somewhat less 
structure on the largest scales, its cosmological parameters 
and other attributes are identical, so our statistical results 
will be similar in both simulations. 

Fig. [8] shows the average fraction of stellar mass in the 



Figure 8. The fractional contributions of instabilities (blue), mi- 
nor mergers (green) and major mergers (red) to the stellar content 
of spheroids, as a function of total galaxy mass, in the Durham 
(solid lines) and MPA (dashed lines) models. 

bulge component of a galaxy that is present due to each of 
the three spheroid-forming mechanisms, as a function of the 
total stellar mass of the galaxy. These fractions include stars 
that may have formed in the bulge as a result of a starburst 
induced by any of the three processes. As expected from our 
previous results, instabilities contribute considerably more 
to bulges in the Durham model, dominating in all but the 
most massive galaxies, where most stars are brought in by, 
or formed in, major mergers. In the MPA model, it is minor 
mergers that are responsible for the bulk of stellar mass in 
spheroids, but, again, major mergers are dominant at the 
highest masses. 

Other galaxy formation models have used slightly differ- 
ent values for and definitions of the mass ratio that ide ntifies 
a major merger. For example ISomerville et all ((2008) use a 
value of 0.25 and include all dark matter within a radius 
of twice the NFW scale-length of the halo. We have tested 
the sensitivity of the above result to the baryonic mass ratio 
used to define a major merger in the Durham model. Val- 
ues between 0.25 and 0.35 modify the fractions in Fig. [S]by 
only a few percent and a value of < 0.01 is required in or- 
der to make major mergers the dominant channel for bulge 
growth above a stellar mass of ~ 5 x 10 10 /i _1 M s . Consider- 
ing instead a versions of the Durham model which employs 
a more realistic prescription for gas stripping in satellites 
l|Font et al.ll200ct ). we find that major mergers occur more 
frequently. This appears to be mainly the result of a higher 
satellite mass (in cold baryons) at the point of merger, which 
reduces the mass ratio. Also, since the total baryonic satel- 
lite mass is larger, the merger timescale is shorter, which 
in turn means that the central galaxy has less time to in- 
crease the mass ratio by cooling more gas before the merger 
occurs. Nonetheless, in this model, major mergers still ac- 
count for < 20% of bulge stars in galaxies less massive than 
~ 5 x 10 11 hT 1 Mq and at higher masses they contribute al- 
most identical fractions to those seen in the Bower et al. 
model. 
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Several authors have considered the role of mergers in 
spheroid formation, although none to our knowledge has 
included disk inst abilities alo n gside mergers as a channel 
for bulge growth. iKang et al.l (|2007l ) used a semi-analytic 
model similar in structure to the Durham and MPA mod- 
els and were able to match the observed fraction of boxy 
early-type galaxies as a function of luminosity and halo 
mass by assum ing that gas poo r majo r mergers result in 
boxy remnants. iKhochfar fcS ilk (2006) used another semi- 
analytic model, though without AGN feedback or N-body 
dark matter merger trees, to examine the contribution of 
'quiescent 'and 'merger' processes - stars formed in progen- 
itor disks and in merger-induced starbursts respectively. 
They found that bulges and elliptical galaxies are dominated 
by quiescently form ed stars. The aspects of sph eroid forma- 
tion in vestigated by IKhochfar fc Silkl |2006h and lKang et al.l 
|2007ft and their methods of analysis differ considerably 
from those presented here, making it difficult to determine 
whether their results are consistent or not with spheroid 
forma tion in the Durham and MPA models. iHopkins et al.l 
(2008) used a 'semi-empirical' model, based on an observa- 
tionally constrained halo occupation distribution, combined 
with theoretical estimates for merger timescales to track 
galaxy mergers. They found that gas rich major mergers 
can account for the observed mass-density of all red, early 
types at z — and the dichotomy in the kinematic proper- 
ties of ellipticals. These results appear, qualitatively at least, 
to be inconsistent with the Durham and MPA models, sug- 
gesting that major mergers alone are sufficient to explain 
the formation of spheroids. The frequency of major merg- 
ers implied by their results certainly seems to exceed that 
in either the Durham or MPA semi-analytic prescriptions, 
though we note that they do not attempt to follow the pro- 
cesses of cooling, star formation and feedback that influence 
the major merger rates in these two models. 

As the above results have demonstrated, not only are 
major mergers a relatively infrequent occurrence in both the 
Durham and MPA models, but their overall fractional con- 
tribution to building spheroids is typically small compared 
to the other two mechanisms in all but the largest galax- 
ies. The importance of disk instabilities, particularly in the 
Durham model, suggests a possible observational test of this, 
perhaps unexpected, conclusion, because the rotation prop- 
erties of spheroids should differ depending on whether they 
formed predominantly through mergers or instabilities. 

The relative importance of rotational and random mo- 
tions in spheroidal systems is measured by the ratio of 
rotational speed to velocity dispersion, V/a, which is of- 
ten plotted against the galaxy's ellipticity to determine 
whether its sh ape is supported by rotation or by veloc- 
ity a n isotropy (llllingworthlll977l; iBrnnevl Il978l; iKormendvj 
19931 ; iKormendv fc Kennicuttl |2004| ; iKormend v fc Fisherl 



2008). During major mergers, tidal forces transfer angu- 
lar mome ntum from the ba r yonic component t o the dark 
halo (e.g.. iFrenk et al.lll985l ; iZavala et ail 120081 ') so the re- 
sulting merger remnants tend to have low V/a and a 
shape that is supported by velocity anisotropy. Spheroids 
forming through disk instabilities on the other hand are 
expected to retain a significant proportion of the angu- 
lar momentum of the progenitor disk, ending up close to 
the 'isotropic rotator' line in the V/a vs flattening dia- 
gram. Our findings are qualitatively consistent with ob- 



servations: low luminosity elliptic als jDavies et alj 19831) 
and m ost 'classical' spiral bulg e s Jlllingworth fc Schechter] 
1 19821 ; IKormendv fc Illingworthl 1 1982ft are found to be 
nearly isotropic, oblate rotators, whereas most giant ellipti- 
cals have insigni ficant rotation (|Bertola fc Capacciolil Il975l ; 
llllingworthlll977M . A detailed analysis of fast and slow rota - 
tors using integral field spectroscopy (|Emsellem et alj|2007f ) 
leads to similar conclusions, with massive ellipticals exhibit- 
ing negligible rotation and the fastest rotators tending to be 
of lower luminosity. 

An important caveat to the test described above is 
that major mergers can also produce spheroids that are 
disky and rapidly rotating. Several simulations of mergers 
between equal mass, gas rich disks have shown that, de- 
spite a large burst taking place, enough gas may survive to 
form a n ew, rotationally supported stellar disk within the 
remnant dBarnes fc Hernquistl Il996l; iBournaud et al.l |20Q. 



Springel fc Hernquist 2005; N aab et al.ll2006l ; IHopkins et al 



2008). Me rger geometry is also like l y to influence remnan t 



structure l|Hernquist fc Mihos|[l995l ; [Bournaud et alj|2005ft . 
These considerations may make our proposed test of the 
models less clear cut. 

The differences in the fractions of spheroid stars added 
by instabilities in the two models are unsurprising given the 
contrasting prescriptions for instability already discussed, 
but it is not simply the dominance of this mechanism that 
decreases the fractional contribution from minor mergers in 
the Durham case. The average mass added by minor mergers 
is typically somewhat larger in the MPA model for any given 
galaxy. This difference appears to stem from two aspects of 
the models. The Durham model has a longer characteris- 
tic star formation timescale and also stronger supernovae 
feedback, which inhibits star formation in lower mass galax- 
ies and efficiently depletes the cold gas available to form 
stars in satellites. Hence, there is less mass on average in 
each Durham satellite. At the high mass end, however, su- 
pernova feedback is less effective, and so the mass brought 
in by major mergers is more similar in the two models. As 
Fig. [5] shows, the two models, in fact, agree rather well on 
the fractional contribution of major mergers to spheroid for- 
mation. 



6 CONCLUSIONS 

We have investigated the origin and evolution of galaxy 
morphology in the ACDM cosmology using two differ- 
ent and publicly available semi-analyti c galaxy formation 
models, the Durham and MPA models l|Bower et al.l 120061 ; 
|Pe Lucia fc Blaizod 120071 ) . both based on the Millennium 
dark matter simulation. Morphology was defined for each 
galaxy by means of the bulge-to-total ratio in K-band lu- 
minosity; we distinguish between the behaviour of bright 
(M K - 51ogh < -22.17) and faint (M K - 51ogh > -22.17) 
galaxies. 

There are many similarities, but also fundamental dis- 
agreements in the predictions of the two models. One of the 
largest differences involves the redshift evolution of the mor- 
phological fractions: spirals are the most common systems 
at all times in the MPA model, while ellipticals dominate 
in the Durham model at redshifts z > 4. The main cause 
of this and other differences highlighted below is the differ- 



Galaxy morphology in the ACDM cosmology 13 



ent treatment of disk instabilities in the two models. In the 
Durham model, instabilities are not only more common at 
early times, but also have more catastrophic consquences: 
the sudden conversion of the entire disk, first into a bar and 
then into a spheroid. Nevertheless, the two models predict 
a fairly comparable morphological mix at the present day. 

In comparison with data from [Benson et all (|2007r ), the 
substantial systematic uncertainties involved in the process 
of fitting bulge and disk light profiles mean that only weak 
constraints can be placed on the fraction of each morphologi- 
cal type as a function of absolute r-band magnitude. Neither 
model is obviously inconsistent with the data, but clearly 
smaller systematic uncertainties would be desirable. The 
bright galaxy populations consist of 51%, 15% and 34% of 
spirals, SOs and ellipticals respectively in the Durham model, 
and 67%, 8% and 25% respectively in the MPA model. 

For typical Milky- Way type haloes, we find that the 
formation time of the halo has, at most, only a small effect 
on the morphology of the central galaxy. In both models, 
haloes hosting ellipticals form slightly earlier than haloes 
hosting spirals or SOs, although this effect is only marginal 
in the Durham model. 

Our analysis confirm s previous studies (|Kauffmannl 
ll996j ; lDe Lucia et al.ll2006h that show that ellipticals galax- 
ies form most of their stars before the bulk of the galaxy is 
assembled. In fact, the ellipticals form and assemble their 
stars earlier than galaxies of other present-day morphologi- 
cal types, although the late 'assembly' times of spirals are a 
consequence of in-situ star formation, rather than mergers. 
In general, the MPA model predicts slightly earlier star for- 
mation for galaxies of all types. In the Durham model, there 
is a significant population of SOs and ellipticals that formed 
at relatively low redshift from the recent collapse of unstable 
disks, yielding spheroids that never regrow a sizeable stel- 
lar disk. The differences in the stellar formation histories 
predicted by the two models are, in principle, accessible to 
observational tests. 

Perhaps the most surprising conclusion of this study 
is that, with the exception of the brightest ellipticals (stel- 
lar mass M* > 2.5 x 10 11 1i~ 1 Mq), major mergers are not 
the primary mechanism by which most spheroids (ellipticals 
and spiral bulges) assemble their mass. Elliptical galaxies are 
more likely than spirals or SOs to have undergone a recent 
major merger, but nonetheless more than half of them do 
not experience one throughout their whole formation his- 
tory (involving galaxies above the resolution limits of our 
calculations). Major mergers are even rarer for spirals, af- 
fecting < 4% of galaxies; their bulges are almost universally 
formed by minor mergers or disk instabilities. These conclu- 
sions would appear to be quite robust for galaxy formation in 
a ACDM universe since they are common to both the MPA 
and Durham models, in spite of the substantial differences 
in the detailed modelling of key physical processes. 

An important difference between the two models con- 
cerns the channel through which less massive ellipticals and 
bulges acquire their mass. In the Durham model, disk in- 
stabilities are responsible for the bulk of the stellar mass in 
spheroids at z = 0, whereas in the MPA model, minor merg- 
ers contribute the most. In the Durham model at least, these 
results appear to be relatively insensitive to the mass ratio 
at which a merger is defined as 'major'. Varying the value 
between 0.25 and 0.35 typically alters the fractions in fig. [8] 



by less than a few percent for any given stellar mass; in order 
for major mergers to become the dominant channel for bulge 
growth in galaxies larger than ~ 5 x 10 10 h~ 1 Mo, the ratio 
would need to be lowered to < 0.01. Adding a more realistic 
prescription for gas stripping of satellites l|Font et alj|2008h 
increases the importance of major mergers, but they still 
contribute < 20% of spheroid mass in galaxies smaller than 
~ 5 x 10 11 /i _1 Mq. Mergers as a whole contribute less mass 
to spheroids in the Durham model because strong supernova 
feedback suppresses star formation especially strongly in low 
mass satellite galaxies. 

The dominant role of disk instabilities in the production 
of spheroids, particularly in the Durham model, offers a pos- 
sible observational test, since spheroids formed this way will 
tend to rotate faster than spheroids formed through merg- 
ers. However, this test would be complicated by the fact that 
gas rich major mergers can also produce high V/a remnants. 
What is clear is that a more precise comparison between the- 
ory and observations will require more detailed modelling of 
disk instabilities than we have attempted here. 

More generally, semi-analytic techniques offer, at 
present, only a somewhat crude way to model the processes 
that establish galaxy morphology. It is therefore remarkable 
that such a simplified approach can lead to some important 
conclusions that appear to be robust to the modelling de- 
tails. This approach also serves to identify specific processes 
that deserve further study, for example, through direct nu- 
merical simulation. 
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